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We review on a chiral power counting scheme for in-medium chiral perturbation theory with 
nucleons and pions as degrees of freedom [[[J. It allows for a systematic expansion taking into 
account local as well as pion-mediated inter-nucleon interactions. Based on this power count- 
ing, one can identify classes of non-perturbative diagrams that require a resummation. We then 
calculate the nuclear matter energy density for the symmetric and purely neutron matter cases up- 
to-and-including next-to-leading order (NLO), in good agreement with sophisticated many-body 
calculations. Next, the neutron matter equation of state is applied to calculate the upper limit for 
neutron stars, with an upper bound around 2.3 solar masses, large enough to accommodate the 
most massive neutron star observed until now. We also apply our equation state to constraint Gn 
in exceptionally large gravitational fields. 



Sixth International Conference on Quarks and Nuclear Physics, 
April 16-20, 2012 

Ecole Polytechnique, Palaiseau, Paris 



* Speaker. 

' I would like to thank my collaborators A. Dobado, A. Lacour, F. Llanes and U.-G. MeiBner for their collaboration 
in several parts of the results presented. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Chiral EFT for nuclear matter 



J. A. Oiler 



1. Introduction 

An interesting achievement in nuclear physics would be the calculation of atomic nuclei and 
nuclear matter properties from microscopic inter-nucleon forces in a systematic and controlled way. 
This is a non-perturbative problem involving the strong interactions. In the last decades, Effective 
Field Theory (EFT) has proven to be an indispensable tool to accomplish such aim. In this work 
we employ Chiral Perturbation Theory (CHPT) to nuclear systems |§, ||[ with nucleons and 
pions as the pertinent degrees of freedom. For the lightest nuclear systems with two, three and four 
nucleons, it has been successfully applied [|5j]. For heavier nuclei one common procedure is to em- 
ploy the chiral nucleon-nucleon potential with standard many-body methods, sometimes supplied 
with renormalization group techniques [^J, or in lattice calculations [p. We have recently derived 
[|l|] a chiral power counting in nuclear matter that takes into account local multi-nucleon interac- 
tions simultaneously to pion-nucleon interactions. Many present applications of CHPT to nuclei 
and nuclear matter only consider meson-baryon chiral Lagrangians (see e.g. [|5|] for a summary), 
without constraints from free nucleon-nucleon scattering. Our novel power counting was applied 
in ref.[§], among other problems, to determine the nuclear matter energy per baryon. We elaborate 
on these results here. Next, we apply them to pure neutron matter case and study the upper limit of 
a neutron star mass. We see that our results can accommodate the recently observed neutron star 
with a mass (1 .97 ± 0.04) M [^J, the largest one confirmed until now. We also use our equation of 
state to constraint the running of the gravitational constant Gn with the gravitational field intensity, 
exceptionally large (around 2 x 10 12 m/s 2 ) inside a neutron star of 2 solar masses. In this case we 



show that Gn cannot exceed its value on Earth by more than a 12% [|10||. 
2. Chiral Power Counting 



Ref. Ql lp establishes the concept of an "in-medium generalized vertex" (IGV). Such type of 
vertices result because one can connect several bilinear vacuum vertices through the exchange 
of baryon propagators with the flow through the loop of one unit of baryon number, contributed 
by the nucleon Fermi seas. At least one is needed because otherwise we would have a vacuum 
closed nucleon loop that in a low energy EFT is buried in the chiral higher order counterterms. It 
was also stressed in ref.[|l2|] that within a nuclear environment a nucleon propagator could have a 
"standard" or "non-standard" chiral counting. To see this note that a soft momentum Q ~ p, related 
to pions or external sources can be associated to any of the vertices. Denoting by k the on-shell 
four-momenta associated with one Fermi sea insertion in the IGV, the four-momentum running 
through the / nucleon propagator can be written as pj = k + Qj. If Q® = &{m n ) = &{p) one 
has the standard counting so that the baryon propagator scales as ff{p^ 1 ). However, if is of 
the order of a kinetic nucleon energy in the nuclear medium then the nucleon propagator should 
be counted as &{jp~ 2 ). This is referred as the "non-standard" case [12]. In order to treat chiral 



Lagrangians with an arbitrary number of baryon fields (bilinear, quartic, etc) ref.[P considered 



firstly bilinear vertices like in refs.[ll , 12], but now the additional exchanges of heavy meson fields 
of any type are allowed. The latter should be considered as merely auxiliary fields that allow one to 
find a tractable representation of the multi-nucleon interactions that result when the masses of the 
heavy mesons tend to infinity. These heavy meson fields are denoted in the following by H, and a 
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heavy meson propagator is counted as due to their large masses. On the other hand, ref.[[l|] 

takes the non-standard counting case from the start and any nucleon propagator is considered as 
&(p~ 2 ). In this way, no diagram whose chiral order is actually lower than expected if the nucleon 
propagators were counted assuming the standard rules is lost. In the following m n ~ k F ~ &{p) 
are taken of the same chiral order, and are considered much smaller than a hadronic scale A x of 
several hundreds of MeV that results by integrating out all other particle types, including nucleons 
with larger three-momentum, heavy mesons and nucleon isobars [Q]. The final formula obtained in 
ref . [[I]] for the chiral order v of a given diagram is 

v* v 

v =4_ £ + £( n . + 4_ 4 ) + £( rf . + v . + fi) ._2)+y p . (2.1) 

(=1 ;=1 

where E is the number of external pion lines, n\ is the number of pion lines attached to a vertex 
without baryons, l\ is the chiral order of the latter with V x its total number. In addition, d\ is 
the chiral order of the i vertex bilinear in the baryonic fields, V, is the number of mesonic lines 
attached to it, ft), that of only the heavy lines, V is the total number of bilinear vertices and V p 



is the number of IGVs. It is important to stress that v given in eq.( |2.1[ ) is bounded from below 



[jjj]. Because of the last term in eq.(2T) adding a new IGV to a connected diagram increases the 



counting at least by one unit. The number v given in eq.(gjj) represents a lower bound for the 
actual chiral power of a diagram, }i, so that jJ. > V. The real chiral order of a diagram might be 
different from v because the nucleon propagators are counted always as 0(p~ 2 ) in eq.(2T), while 



for some diagrams there could be propagators that follow the standard counting. Eq.( |2.1| ) implies 
the following conditions for augmenting the number of lines in a diagram without increasing the 
chiral power by adding i) pionic lines attached to mesonic vertices, l\ = rij = 2, ii) pionic lines 
attached to meson-baryon vertices, dj = v; = 1 and iii) heavy mesonic lines attached to bilinear 
vertices, di = 0, ft), = 1. 

3. Nuclear matter energy density 



The energy per baryon in nuclear matter at NLO in the counting of Eq. (|2JJ) requires to evalu- 
ate the set of diagrams shown in Fig. |l[ Diagram 1 represents the kinetic energy from the Fermi sea 
of nucleons. Diagram 2 arises from the nucleon self-energy summed for all the nucleons. Finally, 
diagrams 3 correspond to the AW interactions in the nuclear medium, 3.1 is for the direct part and 
3.2 for the crossed one. 

A detailed derivation of the final result is given in Ref. [^J. We reproduce here the final 
expression for the diagrams 3, £■$, given by 

^3=4 £ £ (27+ l) X (S£I) 2 f f« d(U - \a + q\)6(U - |a-q|) 

i,j,e,sai,a2 J \ zn ) \ Ln ) 

r Three c 2-2 2n,~ v f d?P ( -|« + Pl)-0(5q2-|»-Pl) y 
x -Tjj(£,£,S;q ,a ,q ) + g Z^ - m J | p2 _ q2 _ /£ Lpi ~~j?\. 

(3.1) 

where an expansion in NN partial waves in the nuclear medium is used. For further details in the 
notation we refer to Ref. [pj. The only quantity not determined from the NN scattering in vacuum 
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Vp = l 

Leading Order 



Next-to-Leading Order 
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Figure 1: Contributions to the nuclear matter energy up to NLO or ^(p 6 )- The wiggly lines correspond to 
the one-pion exchange plus the AW contact interactions from the 0(p°) quartic nucleon Lagrangian [p 





p(fm 3 ) p(fra 3 ) 

Figure 2: Energy per baryon for symmetric nuclear matter (left) and pure neutron matter (right). 



is the constant go. However, one can determine its natural size from the way it is introduced [g], 
around —mm n /4n ~ 0.5 \0m\. The resulting curves for the energy per baryon are shown in Fig. ||, 
left panel for symmetric nuclear matter and right one for purely neutron matter. The preferred final 
values of go are around — 1 m 2 for the former and —0.5 rrv\ for the latter. We see a good comparison 



with sophisticated many body calculations [13] shown by the dotted lines. The experimental point 
for the symmetric nuclear matter panel is indicated by the crossed. In this case we also are able 
to reproduce perfectly the nuclear matter incompressibility K = 259 MeV, to be compared with 
experiment K = 250 ± 25 MeV Q]. 



4. Application to neutron stars. Constraining G 



We want to address two question by employing the energy per baryon (neutron matter equation 
of state at T = 0) obtained in the previous section. The first question is to know whether this equa- 
tion of state, calculated from first principles, is able to account for the large mass of 1 .97(4) M Q of 
the recently observed pulsar J1614-2230 The second question is to know whether the Newto- 
nian constant Gn has the same value as in the Earth for extremely large gravitational fields like those 
found in the previous pulsar, with its estimated acceleration on the surface around 2 x 10 12 m/s 2 
[10]. These questions were addressed in Ref. [[HJ and it was found that our equation of state 
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for neutron matter is able to accommodate such massive neutron stars, with an upper bound for a 
neutron star mass of around 2.3 M Q . This is shown in the left panel of Fig. ||. 

Regarding Gn the procedure was to use the equation of state of Ref. [||] up to Fermi momenta 
of 450 — 600 MeV (leaving this interval as a source of error) and then employ the hardest possible 
equation of state (with a sound velocity equal to the speed of light). The idea is that if the equation 
of state is harder this makes Gn to increase in order to provide enough gravitational attraction. Of 
course, Gn cannot grow indefinitely because otherwise the upper limit value for the mass of the 
neutron star would decrease too much, being in disagreement with the experimental determination 
for the pulsar J1614-2230 [0]. Taking into account all these considerations for the intense gravi- 
tational field in such neutron star, G^ cannot exceed 12% of its value on Earth at 95% confidence 
level. In the right panel of Fig. [5] we show this determination and others from other sources. 





Figure 3: (Color online) Left panel: The neutron star mass as a function of the radius. Right panel: Newto- 
nian constant normalized by its accepted value 6.6738(8)A^(m/%) 2 . From left to right: laboratory on Earth; 
orbital determinations of binary pulsars; white dwarf structure; neutron stars with 1.4 solar masses; neutton 
star with 1.97(4) solar masses. 
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